Energy dependence of nucleus-nucleus potential close to the Coulomb barrier 
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The nucleus-nucleus interaction potentials in heavy-ion fusion reactions are extracted from the mi- 
croscopic time-dependent Hartree-Fock theory for mass symmetric reactions 16 0+ 16 0, 40 Ca+ 40 Ca, 
48 Ca+ 48 Ca and mass asymmetric reactions 16 O+ 40 ' 48 Ca, 40 Ca+ 48 Ca, 16 O+ 208 Pb, 40 Ca+ 90 Zr. 
When the center-of-mass energy is much higher than the Coulomb barrier energy, potentials deduced 
with the microscopic theory identify with the frozen density approximation. As the center-of-mass 
energy decreases and approaches the Coulomb barrier, potentials become energy dependent. This 
dependence signs dynamical reorganization of internal degrees of freedom and leads to a reduction 
of the " apparent" barrier felt by the two nuclei during fusion of the order of 2 — 3% compared to the 
frozen density case. Several examples illustrate that the potential landscape changes rapidly when 
the center-of-mass energy is in the vicinity of the Coulomb barrier energy. The energy dependence 
is expected to have a significant role on fusion around the Coulomb barrier. 

PACS numbers: 25.70. Jj, 21.60. Jz 



I. INTRODUCTION 

Heavy-ion fusion reactions give important informa- 
tion on dynamical evolution and dissipative phenomena 
in a quantum many-body system. Macroscopic mod- 
els [l|, H, H, |3, HI using suitable estimates of nucleus- 
nucleus potentials [3, 0, 0J1, El floL [ill E2I ] , and then 
coupled-channels theories [131 \l4 . 1 1 51 ] have been widely 
used to describe the entrance channel of fusion reac- 
tions. These models underline that the interplay be- 
tween nuclear structure and dynamical effects is crucial 
to properly describe fusion reactions at energies close to 
the Coulomb barrier. While in general rather successful, 
these methods have in common several drawbacks. First, 
nuclear structure and dynamical effects should be treated 
in a unified framework. Second, important effects should 
be guessed a priori. This has been illustrated recently to 
understand new high precision measurements at extreme 
sub-barrier energies [S HE El 0J, HI, where different 
effects like incompressibility |21j . nucleon exchange [22| . 
and the transition from di-nuclear to compound nucleus 
descriptions (23| have been invoked to understand ex- 
perimental observations [H HE El El, HI • Then, the 
hypothesis could only be checked a posteriori. From this 
point of view, it is first desired to use theories where both 
nuclear structure and nuclear dynamics are considered in 
a unified framework. Second, the use of theories where 
all the physical effects mentioned above are automatically 
incorporated can be of particular interest to disentangle 
different contributions. 

Mean-field theories based on Skyrme Energy Density 
Functional (EDF) provide a rather unique tool to de- 
scribe nuclear structure and nuclear reactions over the 
whole nuclear chart. In nuclear reactions, application 
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of the so-called Time-Dependent Hartree-Fock (TDHF), 
more than 30 years ago [H, HI HI HE HI HI HJ, 
to heavy-ion fusion reactions was a major step. With 
the increasing computer power, more and more accurate 
description of the nuclear reactions has been achieved. 
Most recent TDHF simulations include the spin-orbit 
force [HI [H, [H, H HI, HI] and break all the symme- 
tries (plane and axis symmetries generally assumed to 
speed up calculations). Moreover, all the terms of the 
Skyrme EDF used in nuclear structure can now be in- 
cluded [H [H, [M HI HO, E3, El, 13. The possibility 
to perform full three-dimensional calculations and to use 
effective forces consistent with nuclear structure is cru- 
cial to account for the richness of nuclear shapes that 
are accessed dynamically. In addition, the great interest 
of dynamical mean-field theories with respect to other 
methods is that many effects which are known to affect 
fusion such as dynamical deformation, nucleon exchange, 
and nuclear incompressibility are automatically incorpo- 
rated. 

The original purpose of the present work was to bench- 
mark the method proposed in Refs. [H, [33| to obtain 
nucleus-nucleus potential and one-body dissipation from 
the microscopic TDHF dynamics. The possibility to 
obtain such a potential from a mean-field theory has 
been studied in several works using the static EDF tech- 
nique [H, H| HE Ell- More recently, a method called 
density-constrained TDHF (DC-TDHF) [H, which com- 
bines TDHF dynamics with minimization technique un- 
der constraints on the one-body density, has been applied 
in Refs. [UH^], the latter being able to incorporate pos- 
sible dynamical effects through the use of realistic density 
profiles obtained during the evolution. 

Here, we consider a different approach based on 
a macroscopic reduction of the mean-field dynam- 
ics, called hereafter Dissipative-Dynamics TDHF (DD- 
TDHF). This technique could a priori give access not 
only to nucleus-nucleus potential but also to friction co- 
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efficients which play an important role in macroscopic 
models [5(1 [EH H3, [EH and has rarely been obtained from 
fully microscopic theories [1^, [H3| • The main difficulty of 
the macroscopic reduction is to guess the relevant col- 
lective degrees of freedom and their equation of motion. 
Most models assume that the fusion problem can be re- 
duced to a one-dimensional problem on the relative dis- 
tance between nuclei. Here, we will test this hypothesis 
in TDHF and suppose that the dynamics is described 
by a one-dimensional macroscopic dissipative dynamics. 
Since there is a freedom in the choice of macroscopic 
equations, the simple assumption made should first be 
validated. Thanks to alternative techniques used to infer 
Coulomb barriers from TDHF [3(1, [4l| , we show that the 
DD-TDHF method can be a useful tool to get precise in- 
formation on potentials felt during fusion. In particular, 
due to dynamical effects, the deduced nucleus- nucleus 
potentials depend explicitly on the center-of-mass energy 
close to the Coulomb barrier. This energy dependence, 
which was discussed in different models (5f| [56| , is stud- 
ied in detail. 

In this article, we concentrate on nuclear potentials 
study. Aspects related to dissipation will be discussed in 
Ref. (57|. The paper is organized as follows. Next section 
is devoted to the introduction of the DD-TDHF method. 
In Sec. Ill, we give the results and discussions on the 
extracted quantities. A summary is given in Sec. IV. 



II. NUCLEUS-NUCLEUS POTENTIAL FROM 
MICROSCOPIC MEAN-FIELD MODEL 

In this section, the extraction of nucleus-nucleus inter- 
action potentials from mean-field theories is discussed. 
In macroscopic models, these potentials are generally dis- 
played as a function of a few macroscopic collective de- 
grees of freedom describing, e.g., the relative distance, 
shapes of nuclei, mass asymmetry. Here, we will con- 
centrate on head-on collisions between initially spheri- 
cal nuclei 1 and assume that the collective space simply 
identifies with the relative distance R between colliding 
nuclei. The validity of this approximation in the TDHF 
context will be discussed below. In the following, we first 
present a general discussion on different methods to ex- 
tract nucleus- nucleus potential from a mean-field theory. 



A. Some remarks on EDF 

The basic ingredient of a nuclear mean-field model is 
the energy functional of the one-body density p denoted 
by £[p]. In the nuclear context, £[p] is expressed in terms 
of a few parameters generally related to the associated 
effective interaction. Here, we will use the Skyrme EDF 



with the SLy4d [3l| parameters. This choice is partic- 
ularly suited to dynamical calculations because of the 
removal of center-of-mass corrections in the fitting pro- 
cedure of the force parameters [3l|. In the EDF context, 
static properties of nuclei are deduced by minimizing the 
functional with respect to all possible one-body densi- 
ties. The great interest of EDF theory is that the initial 
complicated many-body problem is replaced by an inde- 
pendent particle problem. Indeed, the minimization pro- 
cedure is equivalent to find the set of single-particle states 
that diagonalizes both the self-consistent mean-field de- 
fined through the relation h[p]ij — d£[p]/dpji and the 
one-body density. At the minimum, we have [h[p], p] = 0. 



B. Illustration of fusion with time-dependent EDF 

The static EDF theory has also its dynamical counter- 
part, called time-dependent EDF 2 where the dynamical 
evolution of nuclear systems is replaced by the one-body 
density evolution, i.e., 

m^ = [h[p},p}. (i) 

Dynamical calculations presented in this paper are per- 
formed with the three-dimensional TDHF code developed 
by P. Bonche and coworkers with the SLy4d Skyrme ef- 
fective force [3l[ . As the initial conditions for the TDHF 
time evolution, we prepare the density of colliding nuclei 
by solving static HF equations with the same effective 
force as the one used in the TDHF. The step size in the 
coordinate space is 0.8 fm. Then, we calculate the time 
evolution of the colliding nuclei in the three-dimensional 
mesh. The time step is 0.45 fm/c and the initial dis- 
tance is set between 16 - 22.4 fm, depending on reaction. 
We assume that the colliding nuclei follow the Ruther- 
ford trajectory before they reach the initial distance for 
TDHF. Thus, the initial positions and the momenta of 
the colliding nuclei are determined. As an illustration, 
the density evolution of the 16 O+ 208 Pb head-on collision 
at center-of-mass energy E c . m . = 120 MeV are shown at 
three different relative distances in Fig. [1] 



C. Discussion of nucleus-nucleus potentials 
deduced from mean-field theories 

A good estimate of the interaction potentials felt by 
the two nuclei in the approaching phase within the EDF 
theory could be obtained assuming that the densities of 
the target and projectile remain constant and equal to 
their respective ground state densities. This leads to the 



1 Note that during the reaction nuclei might be deformed. 



2 Although this theory has been called improperly TDHF in the 
past, we will continue to use this acronym in this work. 
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FIG. 1: Top: Density profiles p(X, Y, 0) obtained with TDHF 
for the 16 O+ 208 Pb head-on collision at E c . m . = 120 MeV 
at three different relative distances. The iso-densities (con- 
tour lines) are plotted at each 0.025 fm -3 . The vertical lines 
indicate the positions of separation planes (see text). Bot- 
tom: Total one-dimensional density p(X, 0, 0) (dashed line) 
obtained at the same relative distances. In each case, the two 
solid curves denote respectively pr{X, 0,0) and pp{X, 0,0) 
(see text). Again, the separation plane is presented by the 
vertical line. 



so-called Frozen Density (FD) approximation (see for in- 
stance Ref. In this limit, the interaction potential 
between a target and projectile with their ground state 
densities respectively denoted by pr and pp reads 



1. Definition of the separation plane 

The first step is to properly define the collective coor- 
dinate R that separates the two sub-systems. Here, we 
follow macroscopic models and define the plane of sepa- 
ration at the neck position. In practice, the neck position 
is obtained by considering the two densities 

PT/P (r,t)= J2 l^n /P M)| 2 , 

neT/P 

where ip^{r,t) (resp. ip^{r,t)) denote single-particle 
wave functions initially in the target (resp. projectile) 
which are propagated in the mean-field of the composite 
system up to time t. Note that, if the approaching phase 
is fast enough, these densities have no time to evolve and 
is expected to be very close to their respective ground 
state densities. The separation plane at a given time t 
is then defined at the position where iso-contours of the 
two densities pT(r,t) and pp(r,t) cross. An example of 
densities pp y r(r) (solid lines) as well as the deduced sepa- 
ration plane (vertical thick line) is given for the reaction 
i6 O+ 208p b in bottom of Fig. [T] This figure illustrates 
that the separation plane corresponds to the geometrical 
neck as generally defined in leptodermous systems. 



2. Two-body kinematics 



V* U [R) = £{ Pp+ t}(R) - £[pt] - £\p P ], 



(2) 



where pp+T = Pp + Pt is the total density obtained 
by summing the densities of the target and projectile 
assuming that their centers of mass are at a given relative 
distance R. Note that £[pp+x] here neglects the Pauli 
effect by the two overlapping densities. As we will see the 
following, we use the property of the extracted potentials 
with FD technique essentially at the Coulomb barrier. 
At this point, overlap of the two densities is small and 
Pauli effect is expected to be accordingly small. It should 
however be kept in mind that a proper account of the 
Pauli principle would lead to an increase of the potential 
for small relative distances. 



D. Matching TDHF with a two-body collision 
problem 



Once the separation plane is defined, all quantities rel- 
ative to the dynamics of the two sub-systems can be cal- 
culated. We associate to each sub-space an index i = 1, 2 
and a density pi(t), which equals the total density in the 
sub-space " i" and cancels out in the opposite side of the 
separation plane. Then, all quantities of interest could 
be computed like the number of nucleons in each side of 
the separation plane, i.e., Ai(t) = Tr (pi(t)). Since we are 
considering head-on collisions along x-axis, the center-of- 
mass coordinate Ri{t) and associated momentum Pi(t) 
express as 

Ri(t) = Tr(xpi(t))/Ai(t), Pi(t) = Tr(p x pi(t)). (3) 

We can also compute the inertial mass of the two sub- 
systems, denoted by m„ from the TDHF evolution using 
Tn>i = Pi/Ri. Once these quantities are obtained, the 
TDHF dynamics can be reduced to the two-body collision 
problem where the relative distance R(t) = R% — R2, 
associated momentum P(t) = (m2-Pi— WiP2)/(mi+m2), 
and reduced mass 



The FD approximation is expected to break down if 
strong reorganization of internal one-body degrees of 
freedom occurs in the approaching phase. In the fol- 
lowing, methods to extract nucleus-nucleus potential di- 
rectly from TDHF without assuming frozen densities are 
discussed. 



171x1712 
mi + TU2 



(4) 



are computed at each time step. Note that, when the two 
nuclei are far from each other, the reduced mass properly 
identifies with its initial value p inl — mAj'Ap/{A'p-\-Ap) 1 
where At and Ap denote respectively the initial target 
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and projectile mass 3 and m is the nucleon mass. It is 
worth mentioning that, while for symmetric reactions 
p remains constant during the collision, for asymmetric 
case p may vary after the contact, i.e., p = (a(R). The 
critical discussion on the possible influence of this varia- 
tion on extracted potentials is presented in Sec. IIII El 



E. Dynamical effects on nucleus-nucleus potentials 

The violation of the FD prescription might at least be 
assigned to two effects: (i) In TDHF, the total one-body 
density p(t) does evolve in time. Therefore, the frozen 
density used in Eq. @ should a priori be replaced by the 
density reached dynamically during the evolution to get 
a more realistic nucleus- nucleus potential from TDHF. 
(ii) A second crucial aspect is that part of the relative 
kinetic energy is transformed into internal excitations. 
This effect is generally treated as a dissipative process in 
macroscopic theories. 

In order to include the effect (i) and extract interac- 
tion potentials which account for possible evolution of 
the density along the TDHF path, the so-called density- 
constrained TDHF (DC-TDHF) technique has been de- 
veloped HH E3, Sl| • In this method, at each time-step, 
p(r, t) is deduced from TDHF. Then, the EDF is mini- 
mized under the constraint that the total density matches 
p(r,t). Denoting the minimized energy by £r>c[p(t)](R), 
the potential is then given by 



V^(R) = £ DC [p{t)]{R) - SIM - £[ P p}. 



(•5) 



The great interest of the DC-TDHF method lies in the 
possibility to access the adiabatic potential accounting 
for realistic density profiles. 



F. Matching TDHF with binary dissipative 
collisions 



where V DD (R) and j(R) denote respectively the nucleus- 
nucleus potential and friction coefficient (here "DD" 
stands for Dissipative Dynamics). The friction coefficient 
'y(R) describes the effect of energy dissipation from the 
macroscopic degrees of freedom to the microscopic ones. 
The great interest of this method is the possibility to 
access interaction potentials which account for possible 
dynamical effects and to get information on dissipative 
process from a fully microscopic theory. In the following, 
we show that this method is a valuable tool. 



1. Discussion on the R-dependent mass 

Before presenting applications, additional remarks on 
the origin of the dp(R)/dR term in Eq. ([6]) are manda- 
tory. To obtain the equation, we implicitly assumed that 
the total energy can be written as a sum of a Hamiltonian 
part and a dissipative part E^ ss , i.e. 



E 



V DD (R) + E dil 



2p(R) 



(7) 



where the reduced mass depends explicitly on R. For 
canonical variables (R, P), the Hamilton part in Eq. (|6|) 
arises from the derivative of the kinetic and V DD part in 

previous equation. Using \ R 2 = — ^ ( 2^.) ) > wc 
finally obtain the second equation in ([5]) . The appearance 
of the i?-dependent reduced mass in dynamical mean- 
field calculations has been reported in Ref. [H, and 
turned out to be weak before the Coulomb barrier. In 
all applications presented below we explicitly neglected 
the second term in the evolution of P. Indeed, including 
it or not or taking directly p{R) — Pini has no effect on 
the barrier height and position presented in this work. 
For the sake of completeness this aspect is illustrated in 

sec. mm 



III. APPLICATION 



An alternative tec hniq ue proposed in Ref. [29( , prelim- 
inary tested in Ref. [331 ] . consists in assuming that the 
time evolution of R and its canonical momentum P obey 
a classical equation of motion including a friction term 
which depends on the velocity R: 



dR 
~dt 
dP 
~dt 



P 

_dV DD 
dR 



1 dp{R) • 2 



2 dR 



R z - 1 (R)R, (6) 



3 The discrepancy on the reduced mass found here very close to 
/i(R)//H n i ~ 1 and the one around 0.67 in Ref. [33| was due 
to time-odd terms of the Skyrme field which was not properly 
computed in the original version of 3D code [3l| . This problem 
is now fixed. 



A. Procedure to obtain V DD (R) and 'y(R) 

In Eq. ([S]), we have three unknown quantities: the re- 
duced mass p(R) entering in the first equation, the fric- 
tion coefficient 7(i?) and the nucleus- nucleus potential 
appearing in the relative momentum evolution. 
The reduced mass can be directly deduced from a sin- 
gle TDHF trajectory. Indeed, both R(t) and P(t) can 
be computed at all time along the mean-field trajectory. 
From R(t), R(t) is simply computed from 



R(t) 



R{t + At) - R(t - At) 
2At 



(8) 



where At corresponds to the numerical time step used in 
TDHF. We finally deduced p(R) from P{t)/R(t). 

To obtain the two remaining unknown quantities j(R) 
and V DD (R), a single TDHF trajectory is not sufficient. 
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We assume that the potential energy and the friction 
parameter are not affected by a slight change in center- 
of-mass energy and consider two head-on collisions with 
energies Ei — E c . m . and En — E c . m . + AE (in practice 
AE/E c . m . ~ f — 2% is used). A couple of canonical 
variables (Ri/ii,Pi/ii) is associated to each trajectory. 
Assuming that Eq. © applies in both cases with the 
same potential and friction, we deduce 



j(R) 



[P 



I\Ri=R 



[P 



U\Ri 



[Ri]r t =r - [Rii]r u =r 



(9) 



Then, using one of the trajectories, we obtain dV DD /dR 
as a function of relative distance. The potential V DD (R) 
is deduced by integration over R using its asymptotic 
Coulomb potential at large relative distances. The 
present method clearly relies on the hypothesis that the 
mean-field dynamics could properly be reduced to a one- 
dimensional macroscopic description. As we will see, this 
potential compares rather well with other techniques val- 
idating a posteriori the macroscopic reduction used in 
this work. Finally, V DD (R) is also expected to contain 
dynamical effects like density evolution. 



B. Illustrative example: 16 0+ 16 




FIG. 2: Comparison of potential energies for the 16 0+ ie O 
reaction obtained from different models. The solid, dashed, 
and filled circles-dotted line correspond to the DD-TDHF, 
DC-TDHF [4l|], and FD potentials, respectively. A zoom on 
the Coulomb barrier region is also shown in the insert. 



The potential V DD (R) obtained with the dissipative 
dynamics reduction method for the 16 0+ 16 reaction at 
E c . m . = 34 MeV (and AE = 1 MeV) is displayed (solid 
line) in Fig H The DC-TDHF (dashed line) potential 
obtained in Ref. [4l[ 4 and the FD potential (filled circles- 



dotted line) are also displayed for comparison. Figure [2] 
shows that potentials extracted from the DD- and DC- 
TDHF methods are very close from each other (almost 
identical) even well inside the Coulomb barrier (up to 
R = 5.3 fm). The fact that the potential deduced from 
our method matches the DC-TDHF result gives confi- 
dence in the specific macroscopic equation (Eq. ([5])) re- 
tained to reduce the microscopic dynamics. In addition, 
both methods are almost identical to the FD description 
(for R > 6.5 fm). This indicates that little reorganization 
of densities occurs in the approaching phase. This is in- 
deed confirmed in Fig.[3]where the TDHF density profiles 
obtained at different relative distances (Left) are directly 
compared to densities used in the FD approximation for 
the same R (Right). At and below the estimated bar- 
rier radius Rb = 8.46 fm, little difference between the 
densities Pp/t(X , 0, 0) (Left-solid line) and the ground 
state densities (Right-solid lines) can be seen. As a con- 
sequence, the Coulomb barrier predicted by TDHF is al- 
most identical to the one obtained in the FD case (the 
difference being less than 0.1 MeV). It is worth men- 
tioning at that point that our method assumes neither 
sudden nor adiabatic approximation. Last, another con- 
clusion that could be drawn from the matching between 
DD-TDHF or DC-TDHF and the FD approximation is 
that Pauli blocking effects which are automatically in- 
corporated in the two former approaches and partially 
neglected in V FD do not seem to play a significant role 
close to the Coulomb barrier in 16 0+ 16 0. 



0.2 



0.0 



X 01 
51 

0.0 



0.1 



TDHF density 



Frozen density 



0.0 



■ R = 8'.45 fm ' 


i i i i i i 


R = 7.51fm 

■a/a 




■ R = 6.49 fm 





-10 



5 -10 -5 
X[fm] 



10 



4 Note, however, that a different set of parameters was used for 



the Skyrme effective interaction and numerical aspects. 
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FIG. 3: Left: Density profiles p(X, 0, 0) (dashed) and p P / T (t) 
(solid lines) obtained with TDHF for the head-on le O+ 16 
collision at E c , m , — 34 MeV at three different times. Each 
value of relative distance R is indicated in each left panel. 
Right: Densities p(X, 0, 0) (dashed lines) obtained at the same 
relative distances within the FD approximation are shown. In 
this case, Pp/t(X, 0, 0) (solid lines) identify with the ground 
state densities of the 16 O nucleus. 
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Figure [5] indicates that dynamical effects affect 
marginally the potential felt by the two partners in the 
approaching phase. In the next section, we will indeed 
see that similar conclusions hold in most cases studied 
when J5 c , m . is well above the Coulomb barrier. This is 
the case presented here where E c , m . is three times more 
than the Coulomb barrier. In this limit, the spatial orga- 
nization of nucleons inside each nucleus is almost frozen 
before the contact. 



2 i i i , i , i , i I 

20 40 60 80 100 

V B FD [MeV] 

FIG. 5: Minimal center-of-mass energy, denoted by E^J^ , 
for which the potential deduced from the DD-TDHF method 
identify with the FD case. This quantity is presented as a 
function of the FD barrier energy. 



C. Systematic study of nucleus-nucleus potential at 
high center-of-mass energy 



Previous study is extended to fusion reactions with 
various combinations of nuclei. Since deformed nuclei 
have orientations with respect to the collision axis, which 
increases macroscopic degrees of freedom to be consid- 
ered, we concentrate on collisions involving spherical nu- 
clei. The DD-TDHF technique is applied to the systems 
40 Ca+ 40 Ca, 48 Ca+ 48 Ca for mass symmetric reactions 
and 16 O+ 40 < 48 Ca, 40 Ca+ 48 Ca, 16 O+ 208 Pb, 40 Ca+ 90 Zr 
for mass asymmetric reactions. 
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case. We have checked that this is indeed the case and 
identifies the minimum energy for which the FD limit 
is reached. Potentials obtained with the DD-TDHF are 
displayed by the solid lines in Fig. [4] and are systemati- 
cally compared with the FD approximation. In Figs.[4^a) 
andUJb), the center-of-mass energy used to perform the 
macroscopic reduction is about two times the Coulomb 
barrier energy. Similarly to the 16 0+ 16 reaction case, 
all the examples presented in Figs. [4ja) and B|b) follow 
closely the FD approximation. Higher center-of-mass en- 
ergies have to be used to reach the FD case in the systems 
presented in Figs. BJc) and [5d). We report in Fig. [5] 
the center-of-mass energy threshold, denoted by E^® , 
which corresponds to the minimal E c , m , for which the 
DD-TDHF method gives the FD results (within 5% in 
general) as a function of the FD barrier. The different 
Coulomb barrier energies deduced from the DD-TDHF 
method applied at high center-of-mass energies (denoted 
by Vg D (high E c . m .)) are reported in Tabled These bar- 
riers are systematically compared with the FD case and 
experimental data taken from Refs. [TBI, [HI, [5!|. Over- 
all, we see that the DD-TDHF method applied at high 
center-of-mass energy gives a qualitative agreement with 
experiments. It is however noticeable that the barrier 
height is systematically higher than the experimental ob- 
servation and that the discrepancy increases as ZpZt 
increases. 



FIG. 4: Potentials extracted from the DD-TDHF (solid 
lines) at high center-of-mass energies compared to V FD (filled 
circles- and triangles-dotted lines). 



In all systems, we do expect that, when center-of-mass 
energy increases, the potential will identify with the FD 



We will see in the following that part of the difference 
observed could be understood in terms of departure from 
the FD limit as the center-of-mass energy approaches the 
Coulomb barrier. Indeed, as the energy decreases, densi- 
ties have more time to reorganize. 
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TABLE I: Energy and radii of the Coulomb barrier extracted from the DD-TDHF method. Here, Vg D (high E c .m.) refers 
to barrier deduced for E c . m , > E^.m. while Vg D (low E c , m ,) corresponds to the lowest Coulomb barrier deduced from TDHF 
using Ec.m. — Vg D . The experimental values taken from Refs. [lH l58l. IBl^ are reported when available. 
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D. Center-of-mass energy dependence of 
nucleus-nucleus potential close to the Coulomb 
barrier 

In the previous examples, we have determined the typ- 
ical center-of-mass energy above which the potential de- 
duced from the DD-TDHF technique corresponds to the 
FD case. Here, we show that the extracted potential en- 
ergy is slightly modified as the center-of-mass energy de- 
creases and approaches the Coulomb barrier. As shown 
below, this energy dependence of the nucleus-nucleus po- 
tentials underlines the role of dynamical effects. 



1. Dynamical reduction of the Coulomb barrier energy 



55 




53 1 1 1 1 1 1 1 1 1 1 1 
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FIG. 7: Barrier energy Vg for the 40 Ca+ 40 Ca reaction 
extracted at different center-of-mass energies. The horizontal 
dashed line indicates the FD reference. 
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FIG. 6: Potential energy for the 40 Ca+ 40 Ca reaction ex- 
tracted at different center-of-mass energies. The FD potential 
is displayed with the filled circles-dotted line. 



To illustrate the center-of-mass energy dependence of 
the potential, Figure [6] presents potentials obtained with 
the DD-TDHF method using several center-of-mass en- 
ergies ranging from E c ra = 55 MeV to 100 MeV for the 
40 Ca+ 40 Ca reaction. Again, in the high energy limit, 



potentials identify with the FD case. In addition, an 
increase of center-of-mass energy from E c . m . = 90 to 
100 MeV leads to identical results indicating the sta- 
bility of the method as the energy increases. In oppo- 
site, as E c . m . decreases, potentials deduced from the DD- 
TDHF deviates from the FD case. As E c . m . approaches 
the Coulomb barrier energy, a small change in E c nl . sig- 
nificantly affects V DD as illustrated by the two energies 
E c , m . = 55 MeV and 57 MeV displayed in Fig. [6] In 
order to quantify this dependence, we have reported in 
Fig. [7] values of the Coulomb barrier, denoted by V§ D , 
deduced from the DD-TDHF method as a function of 
center-of-mass energy. Again, if E c . m . is high, Vg D be- 
comes very close to the FD case. As E c . m . decreases, 
Vg D is more and more reduced compared to V£ D . This 
effect, observed in all the cases considered here and called 
hereafter " dynamical barrier reduction" , is a direct con- 
sequence of reorganization of densities in the approach- 
ing phase. This is clearly illustrated in Fig. [8] where den- 
sity profiles obtained for the 40 Ca+ 40 Ca reaction at three 
center-of-mass energies (E c . m . = 55, 57, and 90 MeV are 
shown from top to bottom respectively) and for specific 
R values. In Fig. [51 only the case of E c m . = 90 MeV 
resembles the FD case. At lower energies, a clear de- 
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viation from the FD density profile is observed. Con- 
sidering -E c .m. = 55 MeV, as the two partners approach 
deformation of the two nuclei takes place. This deforma- 
tion initiates the formation of a neck at larger relative 
distances compared to E c , m , — 90 MeV. The center-of- 
mass energy dependence of potential extracted with the 
DD-TDHF technique reflects the difference in the den- 
sity profiles accessed dynamically during the mean-field 
evolution. Note that similar dependence is a priori also 
expected in the DC-TDHF method [U, 5H which ac- 
counts for the dynamical deformation of the densities. 



R = 10.26 fin 



R = 9.82 fm 



R = 9.52 fm 



£ ° 
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„ 5 
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E cm =57MeV 
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H 1 1 1 1 ' 1 1 H 


H 1 1 ' 1 ' 1 ' f— 


-H ' 1 ' 1 ' 1 ' h— 


H ' 1 ' 1 ' 1 ' H 
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X [fm] X [fm] X [fm] 

FIG. 8: Density profiles obtained in TDHF for different rel- 
ative distances R = 10.26 (left), 9.82 (middle), and 9.52 fm 
(right) for the 40 Ca+ 40 Ca reaction at three different center- 
of-mass energies: E c . m , = 55, 57, and 90 MeV from top to 
bottom. These energies corresponds to those used in Fig. [7] 
to obtain V DD (R). 



In all cases considered in this work, a reduction of the 
"apparent" Coulomb barrier seen by the two nuclei be- 
fore fusion is observed compared to the FD case. This 
reduction could always be assigned to large density de- 
formation close to the barrier. In order to quantify the 
magnitude of the dynamical reduction effect, we have 
systematically extracted the lowest barrier energy. This 
quantity (denoted by Vg D (low -E c .m.)), reported in the 
third column of Table U is obtained when the center-of- 
mass energy used in DD-TDHF equals the correspond- 
ing V£ D . In Fig. [9l the difference between the lowest 
barrier and the FD barrier is displayed as a function of 
V£ D . We see that the difference increases almost lin- 
early with V£ D , i.e., with the initial ZpZx- As men- 
tioned previously, the Coulomb barrier energy obtained 
within the FD approximation generally overestimates the 
Coulomb barrier energy deduced from experiments (see 
Table |T|. The discrepancy increases as ZpZx increases. 
Interestingly enough, the lowest energy Vp D is much 
closer to the experimental observation in particular for 
large ZpZ^. 
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FIG. 9: Difference between the barrier obtained with the FD 
approximation and the lowest barrier with DD-TDHF as a 
function of V£ D '. In practice, the lowest barrier is obtained 
by using center-of-mass energy at or close to Vg D . 



2. Critical discussion on the one- dimensional reduction: 
The 16 O+ 208 Pb case 

Previous discussions point out that various density 
profiles might be accessed depending on the center-of- 
mass energy used in TDHF. In macroscopic models, such 
a diversity in densities is usually accounted for by consid- 
ering multidimensional collective space where deforma- 
tion and/or neck are explicitly treated as relevant vari- 
ables [l|, H, H, |3, [H|. Therefore, the energy dependence 
of the potential deduced with the DD-TDHF method 
should a priori be understood as different paths in a 
more complex multidimensional potential energy land- 
scape. As a consequence, one should also a priori con- 
sider macroscopic reduction of TDHF with additional 
collective degrees of freedom which might become ex- 
tremely complicated. 

Here, we show that the simple one-dimensional macro- 
scopic reduction still contains meaningful information on 
the fusion process. We consider the 16 O+ 208 Pb reaction 
for which extensive TDHF calculations have been per- 
formed m m . 

Different potentials deduced for this reaction using 
the DD-TDHF method with different center-of-mass en- 
ergies are displayed in Fig. [TQJ A more complex en- 
ergy dependence is observed in this case compared with 
the 40 Ca+ 40 Ca reaction displayed in Fig. [51 In partic- 
ular at intermediate center-of-mass energies (100 MeV 
< -Ec.m. < 200 MeV), the nucleus- nucleus potential is 
above the FD case. This is clearly illustrated in Fig. [IT] 
where the Coulomb barrier energy V§ D is shown as a 
function of E c _ m ,. A bump for intermediate center-of- 
mass energies is clearly seen. Note that a similar behav- 
ior is also observed in the 40 Ca+ 90 Zr case indicating that 
the energy dependence might be more complex as ZpZx 
increases. However, in view of the potential change com- 
pared to the center-of-mass energy involved (almost two 
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16 Q + 208 pb 



E cm =76MeV 

E cm = 140 MeV 

E am =250MeV 

Frozen density 




extracted. 



FIG. 10: Potential energy obtained with the DD-TDHF 
method for the 16 O+ 208 Pb reaction using different center- 
of-mass energies. The filled circles correspond to the FD ap- 
proximation. 



times Vp D ) this bump is not expected to change dras- 
tically the fusion probability obtained with TDHF. On 
opposite, when the center-of-mass energy is close to the 
Coulomb barrier, a small change in the potential will 
modify significantly the fusion probability. In the follow- 
ing, we will concentrate on this region. 




75 100 



150 200 
E c . m .[MeV] 



250 



FIG. 11: Coulomb barrier energy V B obtained for the 
16 O+ 208 Pb reaction as a function of the center-of-mass en- 
ergy. The horizontal dashed line indicates the FD reference 
while the arrow indicates the low energy fusion TDHF thresh- 
old obtained with TDHF calculations in Ref. 13611. 



The low energy fusion TDHF threshold which is de- 
fined as the minimal center-of-mass energy required to 
fusion in TDHF is also presented as an arrow in Fig. QT] 
A very precise value of this threshold has been ob- 
tained in Ref. 36] using the same TDHF code with 
the SLy4d Skyrme effective interaction by performing a 
large number of TDHF calculations and was found to 
be 74.45 MeV. The lowest barrier energy obtained with 
the DD-TDHF method perfectly matches this threshold. 
Again, this gives additional confidence in this method to 
provide precise information on nucleus-nucleus potential 



E. Effect of coordinate-dependent mass 



In previous sections, the macroscopic equation ^ with 
neglecting the term \^j^-R 2 is used as a starting point 
to extract potentials from the DD-TDHF method. This 
equation is a priori only valid for systems with constant 
reduced mass. This condition is exactly fulfilled by the 
mean-field evolution for symmetric collisions. For asym- 
metric reactions, dependence of the reduced mass with 
relative distance is possible. This situation is illustrated 
in Fig. 1121 where the reduced mass estimated through 
Eq. ([4]) divided by its initial value /ii„i is shown as a func- 
tion of R/R^ D . In all deviation from the initial 
value is observed. In particular, this deviation might be 
significant for the most asymmetric case 16 O+ 208 Pb at 
small relative distances. Similar behavior has been dis- 
cussed in Ref. [13]. In all cases, the center-of-mass energy 
used in the calculation corresponds to V£ D (see table [TJ) . 
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FIG. 12: Reduced mass deduced from the procedure de- 
scribed in Sec . IIII Al divided by its initial value a func- 
tion of relative distance divided by the Coulomb barrier radius 
for the FD case (R^ D ) for the three mass-asymmetric reac- 
tions 16 O+ 40 Ca (a), 40 Ca+ 48 Ca (b) and 16 O+ 208 Pb (c). The 
center-of-mass energies used are reported in each panel. 



To estimate the possible effect of the i?-dependent re- 
duced mass on the extracted potential and friction co- 
efficient, we also extracted potentials and friction coef- 
ficients from Eq. © including the term \ d >1 }^ > R 2 ■ Us- 
ing the same procedure as described in Sec. IIII Ai a new 
potential V DD can be extracted. In Fig. [T3J poten- 
tials with (solid lines) and without (filled triangles) this 
term are compared for the mass-asymmetric 16 0+ Ca, 
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FIG. 13: Comparison of potentials extracted from Eq. © 
(solid line) and from Eq. © without the term ^^^-R 2 (filled 
triangles-doted line) as a function of R/Rg D . The center-of- 
mass energies used are the same as in Fig. [12] 



40 Ca+ 48 Ca, and 16 O+ 208 Pb reactions. In the first two 
cases, the two potentials are indistinguishable while in 
the last case a small difference between the two poten- 
tials is observed. Though we see from Fig. [T^] that the 
reduced mass depends on the relative distance, Fig. [T3l 
shows that this dependence has almost no effect on the 
nucleus-nucleus potential extracted with the DD-TDHF 
method. 



IV. SUMMARY 



automatically incorporate dynamical effects during the 
approaching phase which could be traced back in the en- 
ergy dependence of the nucleus-nucleus potential. This 
energy dependence has been systematically investigated 
for the mass symmetric reactions 16 0+ 16 0, 40 Ca+ 40 Ca, 
48 Ca+ 48 Ca and mass asymmetric systems 16 O+ 40 ' 48 Ca, 

16 O + 208p b; 40 Ca+ 48 Caj 40 Ca+ 90 Zr ^ For tMg system _ 

atic, the following aspects have been discussed: 

• We show that in all reactions the minimal en- 
ergy (Ej[£ ) for which the FD potential is recov- 
ered with the DD-TDHF can always be identified. 
This energy has been systematically investigated. 
We have shown, that E^J^ /V^ D increases as the 
ZpZt increases (Fig. [5]). 

• A clear energy dependence of extracted potential, 
due to dynamical effects which modify density pro- 
file, has been observed in all cases. For systems 
with ZpZt < 400 a continuous decrease of the 
apparent Coulomb barrier is seen as the center-of- 
mass energy decreases (Fig. [7]) while in other sys- 
tems ( 16 O+ 208 Pb and 40 Ca+ 90 Zr) more complex 
energy dependence of the Coulomb barrier energy 
is obtained (Fig. QT]) . 

• In all cases, nucleus-nucleus potential deduced from 
the DD-TDHF method varies rapidly as the center- 
of-mass energy approaches the Coulomb barrier en- 
ergy. Such a rapid change could be assigned to 
the difference in density profiles dynamically ob- 
tained in various TDHF calculations performed 
with slightly different E c/m . 



• Dynamical effects induce a reduction of the appar- 
ent barrier compared to the FD case of the order 
2-3% of V£ D (Fig.©. While the FD Coulomb bar- 
rier generally overestimates the Coulomb barrier es- 
timated experimentally, barriers including the dy- 
namical reduction effect become very close to the 
experimental case (Table QJ. 



The first goal of the present paper was to benchmark 
a technique, called DD-TDHF, to obtain nucleus- nucleus 
potentials and dissipation using a macroscopic reduction 
of mean-field theory based on Eq. ([6]). Several results 
have been obtained that validate the DD-TDHF tech- 
nique: (i) Used with the same conditions, the DD-TDHF 
method leads to potential very close to the DC-TDHF re- 
sult [Il| (Fig. [2]). (ii) As expected, using high center-of- 
mass energies, the DD-TDHF method converges toward 
the FD approximation (Fig. [4j) . (iii) At center-of-mass 
energy close to the Coulomb barrier energy, the dynami- 
cal reduction of the barrier found with the DD-TDHF is 
able to reproduce the low energy TDHF fusion threshold 
obtained in Ref. ^ for 16 O+ 208 Pb (Fig.Qj}. Nucleus- 
nucleus potentials obtained with the DD-TDHF method 



In summary, the DD-TDHF method has been success- 
fully tested in the present work. It gives interesting 
insight in nucleus- nucleus potentials which account for 
dynamical effects. Mean-field calculations show energy 
dependence close to the Coulomb barrier energy. Such 
energy dependence is also expected to affect the sub- 
barrier fusion process. Unfortunately TDHF does not 
provide a correct description of tunneling in collective 
space. The study of sub-barrier fusion with mean-field 
theory would be very interesting but clearly requires to go 
beyond TDHF. Finally, we would like to mention that the 
DD-TDHF technique also gives dissipative kernels from 
a fully microscopic theory. This aspect, which is crucial 
in macroscopic theory, will be discussed in a forthcoming 
article [13]. 
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